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Abstract 

The properties of antiferromagnetic Heisenberg S = \ ladders with 
2, 3, and 4 chains are expanded in the ratio of the intra- and inter- 
chain coupling constants. A simple mapping procedure is introduced 
to relate the 4 and 2-chain ladders which holds down to moderate val- 
ues of the expansion parameters. A second order calculation of the 
spin gap to the lowest triplet excitation in the 2- and 4-chain ladders 
is found to be quite accurate even at the isotropic point where the 
couplings are equal. Similar expansions and mapping procedures are 
presented for the 3-chain ladders which are in the same universality 
class as single chains. 
PACS numbers: 
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I. INTRODUCTION 



The study of antiferromagnetic Heisenberg spin-| ladders has recently become 
of great interest because they offer the possibility of realizing spin liquid states with 
short range RVB (resonance valence bond) character [1-5]. Experimental systems 
currently being investigated are (VO)2P2C>7 [6,7] and the homologous series oi 
compounds Sr„_ 1 Cu rl+1 02n [8,9]. The main evidence to date comes from numerical 
investigations using both Lanczos diagonalization of finite ladders [1,2,10], and 
density-matrix renormalization group [5], and quantum transfer-matrix methods 
[11], These studies have shown that the low energy properties of ladders are 
governed by fixed point behavior determined by the limit of strong interchain and 
weak intrachain coupling. In this paper we present analytic expansions about this 
strong coupling limit for ladders with up to four legs (or chains). We compare these 
expansions with the numerical results and find that they hold qualitatively but 
not always quantitatively down to the isotropic limit (equal inter- and intrachain 
couplings). Further we introduce a mapping scheme to relate the 4-chain ladder 
to the simple (i.e. 2-chain) ladder and the 3-chain ladder to the single chain. 

The Hamiltonian of such a ladder takes the form 

i7 = £jS,-S m + £j ± S,-S m . (1) 
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The intrachain coupling is J and the interchain coupling across the rungs is J±. 
The summation over the length of the chains — which will eventually tend to 
infinity — is denoted by <-> and the summation over the number of chains running 
parallel is denoted by J. 

In the strong coupling limit J± 3> J, the part of the Hamiltonian with the 
sum along the chains is treated as a perturbation of the system of uncoupled 
rungs. For the simplest spin ladder (2-chains) the strong coupling perturbation 
expansion up to third order in J/ Jx is performed in Section2. Ladders with three 
chains coupled together will be discussed in Section 3 and a mapping to the single 
Heisenberg spin-^ chain with an effective coupling determined by the splitting 
between the ground state and the first excited state of a 3 X 2 spin cluster. Lastly 
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the case of 4-chain ladders is treated in Section 4 using two different approaches. 
First the four chains are considered as a system of two simple ladders running 
parallel, with a small coupling J' ± between the two adjacent chains. Second, the 
four chains are treated in the same way as the simple ladder by starting in the limit 
of non-interacting rungs, each with four spins. In the second case — in analogy to 
the 3-chain system — a mapping to the simple ladder is performed. 



II. TWO COUPLED SPIN-i CHAINS: THE SIMPLE SPIN LADDER 

The best studied system is the simple spin-| ladder, i.e. two parallel spin-| 
chains strongly coupled together. The Hamiltonian has the form 

H — H° + H 1 , (2) 

where H° and H 1 are given by 

H° = E K = £ J ± (S\ ■ 8%) (3) 

R=l R=l 

H I =jbnk=jbj(.8R-SM + SR-sUi)- ( 4 ) 

R=l R=l 

The intrachain (interchain) coupling is given by J (respectively Jx) , while 
denotes the spin operator on the R-th rung on chain i. For chains of length L 
periodic boundary conditions are introduced by defining S^ +1 = S\. 

The eigenstates of the rung Hamiltonian are given by a singlet state \s) with 
energy E, — — f Jx and three triplet states \t") with spin z-component a — —1, 0, 1 
and energy E t — |Jx- The eigenstates of H° are direct products of rung states. 



A. Spin Gap and Dispersion relation of magnon excitations 

When applying perturbation theory in J/Jx for the strong-coupling limit, the 
first order correction to the ground state energy vanishes, since H 1 pairwise excites 
two adjacent singlets of the unperturbed ground state |0) = to a linear 

combination of triplets with total spin 5 = 0: 
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h ! \o) = \j E 

Z R 



.{t°t°- t+t~ -t-t+)...s). (5) 



This excitation leads to a correction of the ground state energy. Up to third order 
perturbation in the strong coupling limit, the ground-state energy per rung is 



E *- J 



3 3 

4 + 8 \,.<: 



(j, ) + 16 (j, ) 



(6) 



The first excited state of the unperturbed ladder is obtained by promoting one 
rung to a triplet state. The X-fold degeneracy of this state is lifted in first order 
perturbation in the strong-coupling limit. To first order in J/J± the eigenstates 
are given by Bloch states 

l^*) = 4=E e<fa 'l"---*f ■■■"). ( ? ) 

where the l-th rung is excited to a triplet with S z = a . The excitation energy w(k) 
for magnon excitations up to third order Jj Jx in the strong-coupling limit is 



— ( — I (2 cos k + 2 cos 2k — cos 3k — 3). 
8 V J ± ) 

To second order this result is equal to that obtained by Barnes et. al. in [2], except 
for the additional second order term oc cos2&. The energy has a minimum at 
k — it. Therefore the spin gap A — oj(tt) is 

£=>-(iH(i)'+i(f)* <•> 

In Fig. 1 the dispersion relation in second and third order is plotted for a ratio 
of J/Jx — 0.5. The third order correction improves the agreement to numerical 
data from [2], The deviation near k — is due to 2 magnon processes. Fig. 2 
shows the spin gap in second and third order compared to numerical data for a 
2x8 Heisenberg ladder [1]. At the isotropic point J/Jx — 1 results from Barnes et. 
al. [2] for 2x8 and 2 X oo ladders are included. In the isotropic region, the third 
order correction unfortunately leads to worse agreement than the second order 
correction, which is surprisingly good in the infinite-length limit for the simple 
Heisenberg spin ladder. 



III. LADDERS WITH THREE SPLN-± CHAINS 



A system of three coupled spin-1 chains has a degenerate ground state in the 
strong-coupling limit which leads to an additional complication compared to the 
study of the simple ladder (or any other system of an even number of chains coupled 
together). White et. al. [5] give an explanation of the fundamental difference of 
even and odd number of spin chains coupled together as due to the behavior of 
topological spin defects. Since the rung states are already degenerate, perturbation 
theory for degenerate systems must be used from the beginning. This can be 
realized by mapping the system to the single Heisenberg spin-^ chain — which has 
been studied extensively — using an effective coupling in first order perturbation 
in J I J±. Later we use the exact diagonalization of a 3 X 2 spin cluster to determine 
the effective coupling. 

The Hamiltonian of the system is given in analogy to the one of the ladder by 

H^^ + H 1 , (10) 

with 

tf°=E*£=E ^(sk-si + sj-s^) (ii) 

R-l R-l 

h'= E#* = E £jsj.sj +1 . (12) 

R-l R-l j=l 

The energy levels of a rung system together with their degeneracy are depicted 

in Fig. 3a. The two spin-j states are denoted by \d") (doublet), a corresponding 

to S z = i^, the spin-| state by (quartet) with a corresponding to S z = 

_3 l l 3 

2 > 2 ' 2 ' 2 ' 

The energy levels of an uncoupled pair of rungs (with Hamiltonian H° = 
fljj-! + Hr_ 2 ) and their splitting into spin-subspaces upon switching on the per- 
turbation H 1 — Hj t _ 1 are shown in Fig. 3b. 

In first order perturbation theory the ground state of a pair of rungs is 

E g = -2J ± - | J (13) 
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and the splitting to the lowest triplet is equal to J. The effective coupling in first 
order perturbation theory in the strong coupling limit then is 



Jeff = J (W) 

The results of the exact diagonalization of a 3 X 2 cluster are shown in Fig. 
4. The separation of the second from the first excited state is greater than the 
splitting J e ff of the singlet and the triplet, reaching its minimum of = J e g at the 
isotropic point J/Ji = 1. Thus the mapping should give reasonable results for 
temperatures k&T < J e ff. At isotropy, from Fig. 5 we obtain a value J c g — 0.82J. 
A more complete density matrix renormalization group (DMRG) study of isotropic 
Heisenberg coupled chains by White et. al. [5] leads to qualitatively comparable 
results. They consider finite ladders with open boundary conditions so that there is 
a finite excitation energy determined by the velocity of the des Cloizeaux- Pearson 
mode. In a single chain this is proportional to the coupling constant. The ratio 
of these velocities in the 3-chain and 1-chain systems gives a direct measure oi 
renormalized coupling J e g sa 0.68. 

IV. LADDERS WITH FOUR COUPLED SPIN-i CHAINS 

To study the ladders with 4 parallel coupled Heisenberg spin-| chains in the 
strong coupling limit, two different approaches are taken. First the 4 chains will 
be treated as two simple ladders coupled together and expanded around the strong 
rung coupling limit. Secondly, we map the 4-chain system to a simple double-chain 
ladder with renormalized coupling constants. The Hamiltonian of this system is 
the sum of three terms H = H° + H 1 + H 2 with 

ff° = E J±{S 1 R S^ + S^.S^) (15) 

R=l 

ff 1 = E E J*r-% +1 (i6) 

R=l j = l 

H 2 = E J'± S« ■ S 3 R . (17) 

R-l 
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As in the double chain system, the interchain coupling (coupling between chains 
one and two and between three and four) is J± and the intrachain coupling is J. 
Periodic boundary conditions again are introduced by S* L+1 = S\. The coupling 
between the two ladders is J\ . 



A. Two coupled spin-^ ladders 

The part 7?° + H 1 of H is the sum of two simple ladder Hamiltonians: 



H° 



H 1 = H u 



H,. 



(18) 



The indices u and I denote simple ladder Hamiltonians (2) on the upper resp. lower 
simple ladder. 

To determine the ground state energy and the energy of the low-lying excita- 
tions, perturbation in J and J' ± in the strong coupling limit J J' ± <C J± up to 
second order is applied. The first excited state is 2 ■ 3X-fold degenerate, the factor 
2 arising since promoting a singlet to a triplet on one rung can be done on either 
of the two coupled ladders. The X-fold degeneracy is lifted as in the simple ladder 
system by the perturbation H 1 (transforming to Bloch states); the remaining 3-fold 
degeneracy is a spin degeneracy. The second perturbing Hamiltonian H 2 will lift 
the 2-fold degeneracy into even and odd parity states, the parity-transformation 
-Pchain being defined by reversing the chain order (P c hain : i — > 5 — i). 

For the ground state the energy per rung to second order in J and J' ± is 

2" 



= -J± 



3 3 

5 + 1 



(0 



+ 



32 V J 



(19) 



With the corrections obtained for the first excited state with odd ( — ) and even 
parity (+), the dispersion relation for the 4-chain system as two coupled ladders 
up to second order perturbation in the strong coupling limit is 



uyf{k) 
J± 



1 + 



(£)■ 



f I (3 



,2i)-i(i 

' 32 1 J ± 



(20) 



I (II 
4 I J , 



8 I J, 



i j;j 



cos k 
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The dispersion relation for the two branches is shown in Fig. 6 for J = J' ± = ^J±- 
Fig. 7 shows the spin gap for the two branches. 



B. Four coupled spin-^ chains 

The Hamiltonian is again of the form (15)-(17). The interaction Hamiltonian 
H 2 between ladders will not be treated as perturbation of the one-rung eigenstates 
of the double ladder but instead exact eigenstates of the sum H a -\-H 2 are the basis- 
states for a perturbative treatment of the intrachain coupling J. Additionally the 
problem is mapped to the simple ladder by exact diagonalization of the 4x2 
cluster. 



1. Exact one-rung eigenstates 
The 16 one-rung eigenstates of H° + H 2 are denoted by 





E, T 




Et a 




E ^ 




E h 






«*> 


E q -- 
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There are two singlets \s T ), three triplets \f a ), and with S z components 
a — —1, 0, 1 and a quintet with S z components cr = —2, —1, 0, 1, 2. 
For J' ± — Jx, the ground state energy per spin of the 4-chain system is 



Es 



-Jx 



3 + 2^3 16 + 3^ tJ_ 



24 



VJx 



(22) 



up to the second order in J/J±. By the same reasoning as for the double ladder, 
the coupling leads to a splitting of the first excitation from the ground state 
into an even and odd parity state. For the odd parity state the dispersion relation 



for the low-lying excitations of four coupled Heisenberg spin-| chains up to second 
order perturbation theory in the strong coupling limit is 

U ~ {k) = 0.659 + 1.075 ( J 



J 

J 

h 



(j?-) (1.086 - 0.035 cos k- 0.469 cos 2k), (23) 



OJ+(ifc) /J'. 

— — = 1.366 + 0.667 — cosfe 



J \ 2 

— (1.826 - 1.155 cos k - 0.081 cos 2k) (24) 

J A. J 

for the even parity state. The dispersion relation for the two branches is shown 
in Fig. 8. For the odd parity branch, the spin gap minimum is always at k = 7r. 
On the other hand, the minimum of the even parity branch jumps from k = it to 
k = at J/Ji = 0.6. This explains the kink in the curve of the even parity spin 
gap in Fig. 9. 



C. Mapping to an effective simple ladder 

As in Section 3, let us map the 4-chain ladder into an effective simple ladder to 
determine the spin gap. Throughout this section, we assume J' ± — Ji again. To 
this end, we calculate the eigenenergies of a 4x2 cluster, and determine the effective 
coupling constants of the effective simple ladder, J* — J*( J, Jj_), J ± = J±{Ji J±)> 
so as to reproduce the same low-energy spectrum. In this approximation, the spin 
gap of the 4-chain ladder is then given by that of the simple ladder as 

A4-ch„i„( J J±) = A 2 _ chain ( J*( J, JO, J* (J, JO). (25) 

The energy spectrum of a 2 X 2 cluster is easily calculated. The ground state is 
a spin singlet, and there are two magnons, bonding and antibonding combinations 
of the two states in which one rung is a spin singlet and the other is a triplet. The 
energies of these three lowest states are 
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So = -1(J* + JI) - -J J* 2 - J'J' ± + Jf, (26) 

Ei,± = -±J ± ±±J*. (27) 
These three states are sufficient to determine the coupling constants: 

J* = E 1}+ -E 1} _ , (28) 



Jl = \/D 2 -DJ* -\J* 2 + \r, (29) 

where D = \{E 1A + E t _) - E is the average energy separation between the 
ground state singlet and the two triplets. There is one other spin triplet consisting 
of two magnons, with the energy, E lttt = j( J* — J*). 

Fig. 10 shows the energy levels of the 4x2 cluster as a function of the intrachain 
coupling, J. The ground state is a spin singlet and there are two triplets above 
it. The two triplets are split with a separation increasing with J. This splitting 
corresponds to the band width of propagating magnons along the chain direction. 
At around J/Jx ~ 0.64, another triplet goes down and becomes lower than . 
This triplet is the two magnon state. Therefore, this crossing corresponds to the 
point where the bottom of the two-magnon continuum becomes lower than the top 
of the one-magnon mode. However, these two triplet eigenstates, having different 
parities with respect to the mirror symmetry, -P c h a in, do not mix to each other, and 
this level crossing does not have significant consequences. 

The energy spectrum Fig. 10 has the same structure as the 2x2 system, at 
least for small J's. The corresponding effective couplings are determined by (28) 
and (29). The result is shown in Fig. 11. The mapping breaks down for large 
couplings, J/J± > 0.70, where the condition D > v ^ +1 (E\ : + — Si,-), which is 
necessary for a real J*, is no longer satisfied. 

We can now estimate the spin gap of the 4-chain ladder by using (25). Here 
we use a Pade approximation for the A 2 - C hain determined by the strong coupling 
limit (9) and the weak-coupling asymptotic form, A oc J± [12]: 

A 2 - chflin (J,J ± ) = JiG(^), (30) 

(2 + 2a)z 2 + (1 - 3a)x + 2a 
1 >~ (2 + 2a)x> + (3-a)x + 2 ' 1 > 
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where a = limj L _ J±- This Pade approximant has a correct asymptotic 

form in both limits, J±/J — > oo and J±/J — > 0. We use here a = 1 determined by 
^2-chain — 0.5 at Jx = J, but this form does not agree so well with the numerical 
results [2] at < Ji/J < 1. 

Figure 12 shows the result of the spin gap for the 4-chain ladder. The gap is 
obtained by using (31) for (25) with the values of J*'s and J*'s shown in Fig. 11. 
The values determined by numerical diagonalization [5,10] are included at J± — J. 
The result of the second-order perturbation (23) is also included. This mapping 
shows a correct tendency. 

V. CONCLUSIONS 

In this paper we presented the results of analytic expansions about the limit 
of strong interchain coupling. It is convenient to expand around this strong cou- 
pling limit since the qualitative behavior and therefore the universality class of the 
system remains unchanged. An analytic expansion allows one to examine the evo- 
lution of the system to the case of isotropic coupling. A magnon is an elementary 
S — 1 excitation, but the spin density distribution which is localized on a single 
rung in the strong coupling limit evolves continuously into a distribution spread 
over several rungs at isotropy. With decreasing interchain coupling the magnon 
spectrum changes its dispersion relation from a simple cosine band as the size of 
the magnons increases and longer range hopping matrix elements enter. An 5 = 1 
magnon can be regarded as a bound state of two S = j spinons [13] on individual 
chains, but since the interchain coupling changes the excitation spectrum com- 
pletely this analogy is only qualitative at best. The 4-chain ladder also has a spin 
gap to S = 1 magnon excitations and so belongs to the same universality class as 
the 2-chain ladder, and a mapping procedure between 4-chain and 2-chain ladders 
is possible. This works well down to moderate values of the intra- to interchain 
coupling but not down to the isotropic limit where they are equal. On the other 
hand a second order perturbation for the spin gap works surprisingly well down to 
the isotropic limit. The 3-chain ladder can be mapped onto the single chain and a 
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simple mapping procedure is found to work quite well down to the isotropic limit 
when compared to numerical results. The low energy behavior of the two systems 
will be similar and the spinons will now extend over the three chains. Lastly the 
success of the expansion around the strong coupling limit of the Heisenberg lad- 
ders encourages one to consider a corresponding expansion for the doped ladders 
described by a t-J model which have shown interesting results in mean field and 
numerical investigations [14,15]. 
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FIG. 3. (a) The rung energies for the 3-chain ladder, (b) These combine to give 
the eigenenergies of a pair of rungs, and under the action of H 1 they split yielding the 
spectrum for a 3 X 2 cluster. 




FIG. 4. Diagonalization of 3 X 2 cluster classified by total spin S. 
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FIG. 9. Spin gap for a magnon of four Heisenberg spin-i chains. 
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FIG. 11. Effective couplings when the 4-chain ladder is mapped to the simple ladder. 
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FIG. 12. Spin gap of the 4-chain ladder calculated by mapping to a simple ladder. 
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